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Based on 106 million i/;(3686) events collected with the BESIII detector at the BEPCII, the decay 
V>(3686) — > uKKn is studied. Enhancements around 1.44 GeV/c 2 and the /i(1285) are observed 
in the mass spectrum of KKtt, and the corresponding branching fractions are measured, as well as 
the branching fractions of t/>(3686) — > ujK* + K~ + c.c. and t/j(3686) — > luK*°K° + c.c, all for the 
first time. 

PACS numbers: 13.85.Hd, 25.75Gz 



I. INTRODUCTION 



Charmonium decays play an important role in the study of the strong interactions. The t/>(3686), which was the 
second charmonium state discovered [1], has been extensively studied both experimentally [2-6] and theoretically [7- 
10]. Pcrturbativc QCD [11, 12] predicts that the partial widths for J/ifj and 7^(3686) decays into an exclusive 
hadronic state h are proportional to the squares of the cc wave-function overlap at zero quark separation, which are 
well determined from the leptonic widths. Since the strong coupling constant, a s , is not very different at the J/ip 
and 7^(3686) masses, it is expected that the J/ip and t/;(3686) branching fractions of any exclusive hadronic state h 
are related by 

B(^(3686) -» h) ^ /3(t/,(3686) -> e+e~) ^ 
Vh B{J/i>^h) ~ B(J/^^e+e-) /0 ' 

This relation defines the " 12% rule" , which works reasonably well for many specific decay modes. A large violation 
of this rule was observed by later experiments [3, 4, 6], particularly in the pn decay. Recent reviews [9, 10] of relevant 
theories and experiments conclude that current theoretical explanations are unsatisfactory. Clearly, more experimental 
results are also desirable. 

A pseudoscalar gluonium candidate around 1.44 GeV/c 2 , the so-called £//i(1440), was first observed in pp annihila- 
tion in 1967 [13]. Studies in different decay modes revealed the existence of two resonant structures, the ?y(1405) and 
the 77(1475) [14]. The 77(1475) could be the first radial excitation of 77 (958). The L3 measurements of the KKtt and 
?77r + 7r _ decays in photon-photon fusion suggest that the 77(1405) has a large gluonic content [15]. However, CLEO 
with a five times larger data sample did not confirm the L3 results, but their upper limits are still consistent with the 
glueball and the radial excitation hypotheses for the 77(1405) and 77(1475) [16]. 
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In this paper, the first observation of enhancements at around 1.44 GeV/c 2 and at the /i(1285) resonance in the 
mass spectrum of KKtt (KgK + ir~ +c.c. and K + K~ir°) produced in -0(3686) — > ljKKtt and the measurements of the 
corresponding branching fractions are reported. In addition, the branching fractions of "0(3686) — ^ wK* + K~ -\-c.c. and 
0(3686) — > uiK* K°+c.c. are also measured for the first time. The analysis reported here is based on 1.06 x 1O 8 0(3686) 
events collected with the BESIII detector at BEPCII. 

II. BESIII DETECTOR 

BEPCII/BESIII [17] is a major upgrade of the BESII experiment at the BEPC accelerator [18, 19] for studies of 
hadron spectroscopy, charmonium physics, and r-charm physics [20]. The design peak luminosity of the double-ring 
e+er collider, BEPCII, is 10 33 cm^s" 1 at a beam current of 0.93 A at the ^(3770) peak. The BESIII detector with a 
geometrical acceptance of 93% of 47r, consists of the following main components: 1) a small-celled, helium-based main 
draft chamber (MDC) with 43 layers. The average single wire resolution is 135 fim, and the momentum resolution for 
1 GeV/c charged particles in a 1 T magnetic field is 0.5%; 2) an electromagnetic calorimeter (EMC) made of 6240 Csl 
(Tl) crystals arranged in a cylindrical shape (barrel) plus two endcaps. For 1.0 GeV photons, the energy resolution 
is 2.5% in the barrel and 5% in the endcaps, and the position resolution is 6 mm in the barrel and 9 mm in the 
endcaps; 3) a Time-Of-Flight system (TOF) for particle identification composed of a barrel part made of two layers 
with 88 pieces of 5 cm thick, 2.4 m long plastic scintillators in each layer, and two endcaps with 96 fan-shaped, 5 cm 
thickness, plastic scintillators in each endcap. The time resolution is 80 ps in the barrel, and 110 ps in the endcaps, 
corresponding to a K/tt separation better than a 2 a for momenta below about 1 GeV/c; 4) a muon chamber system 
(MUC) made of 1000 m 2 of Resistive Plate Chambers (RPC) arranged in 9 layers in the barrel and 8 layers in the 
endcaps and incorporated in the return iron of the superconducting magnet. The position resolution is about 2 cm. 

The GEANT4-based simulation software BOOST [21] includes the geometric and material description of the BESIII 
detectors, the detector response and digitization models, as well as the tracking of detector running conditions and 
performance. 1.06 x 10 8 inclusive Monte Carlo (MC) events are used in our background studies. The production of 
the ^(3686) resonance is simulated with the event generator KKMC [22, 23], while the decays are generated with 
EvtGen [24] with known branching fractions [14], and by Lundcharm [25] for unmeasured decays. The analysis is 
performed in the framework of the BESIII Offline Software System (BOSS) [26] which takes care of the detector 
calibration, event reconstruction, and data storage. 
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III. EVENT SELECTION 

In this analysis, the w meson is reconstructed in its dominant decay uj — > 7r + 7r~7r°; Kg and 7r° are reconstructed 
from the decays Kg — > ir + n~ and ir° — > 77. The final states of -0(3686) — > ujK s K + tt~ and u>K + K~ir° arc 
2(7r + 7r _ )A' + 7r _ 77 and TT + n~ K + K ^7777, respectively^ . 

Charged tracks are reconstructed from MDC hits. Each charged track (except those from K s decays) is required to 
originate from within 2 cm in the radial direction and 20 cm along the beam direction of the run-by-run-determined 
interaction point. The tracks must be within the MDC fiducial volume, |cos#| < 0.93, where 9 is the polar angle. 
The information from the TOF and dE/dx is combined to form a probability Prob(A) (Prob(7r) or Prob(p)) under a 
kaon (pion or proton) hypothesis. To identify a track as a kaon, Prob(A) is required to be greater than Prob(7r) and 
Prob(p). 

Electromagnetic showers are reconstructed from clusters of energy deposits in the EMC. The energy deposited in 
the nearby TOF counters is included to improve the reconstruction efficiency and the energy resolution. A photon 
candidate is a shower in the barrel region (| cos#| < 0.8) with an energy larger than 25 MeV, or in the endcap region 
(0.86 < I cos#| < 0.92) larger than 50 MeV, where 9 is the polar angle of the shower. The showers close to the gap 
between the barrel and the endcap are poorly measured and are thus excluded from the analysis. Moreover, the 
EMC timing, with respect to the collision, of the photon candidate must be in coincidence with collision events, i.e., 
< t < 700 ns, to suppress electronic noise and energy deposits unrelated to the event. 

A!j candidates are reconstructed from secondary vertex fits to all oppositely charged track pairs in an event. The 
combination with an invariant mass closest to the nominal K s mass (m K o) is kept. The reconstructed K s is used as 
input for the subsequent kinematic fit. 

IV. 0(3686) cjK^K+tt- + c.c. 

Event candidates should have a K s , four charged tracks with zero net charge, and two or more photons. At least 
one charged track is positively identified as a kaon. When more than one kaon is identified, the one with the maximum 
Prob(A) is chosen. The 0(3686) — > ir + ir~ K s K + tt~^^ candidates are subjected to a four-constraint (4C) kinematic 
fit provided by four-momentum conservation to reduce background and to improve the mass resolution. For events 



1 The charge-conjugate final state is included throughout the paper unless explicitly stated. 
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with more than two photons, the combination with the minimum \ 2 value of the fit (denoted as x\c) * s retained, and 
xlc ^ s required to be less than 40. For ir° candidates formed from a photon pair, the invariant mass of the photon pair 
must be within the range |M 77 — TO^o) < 0.02 GeV/c 2 . The combination of 7r + 7r _ 7r° with an invariant mass closest 
to the nominal u> mass (m u ) is chosen as an u> candidate. Figure 1(a) shows the distribution of the M T +,- versus 
the M„+„-„o invariant mass, where aw — Kg cluster corresponding to "0(3686) — > LuKgK + ir~ is apparent. Events 
are kept for further analysis if the tt + tt~ invariant mass is in the range 0.489 GeV/c 2 < M 7T + 7T - < 0.505 GeV/c 2 and 
the tt + tt~tt° mass in the range 0.743 GeV/c 2 < M^+^-^o < 0.823 GeV/c 2 . Events in the K a s sideband range (0.012 
GeV/c 2 < |M w+7r - - m K o \ < 0.02 GeV/c 2 ) or the w sideband (0.06 GeV/c 2 < \M w +„-„o - m u \ < 0.10 GeV/c 2 ) arc 
used to estimate the background. In addition, to veto the 0(3686) — > ir + TT~J/ip background events, the recoil mass 
against tt + tt~ is required to satisfy M£+ c °0 — mj/^\ > 0.007 GeV/c 2 . 

A. Branching fractions for "0(3686) -¥ cjA(1440) -> ujK s K + n~ + c.c. 

Figure 1(b) shows the invariant mass M K o K +^- for selected events, where a peak around 1.44 GeV/c 2 (denoted 
as X(1440), since we do not distinguish 7/(1405) and ?y(1475) from this analysis) and the /i(1285) are evident. To 
verify that the observed peaks originate from the process "0(3686) — > LuKgK + ir~ , the backgrounds are investigated 
from both data sideband and inclusive MC events. The non-_fT<j and non-w backgrounds, estimated by using events 
in the A"^ and uj sideband regions, and normalized according to the ratio of MC events falling into the sidebands 
to the signal region, are shown as the shaded histogram in Fig. 1(b). No evident peak around 1.44 GeV/c 2 is seen. 
Other potential ^i>(3686) decay backgrounds are checked with 106 million 0(3686) inclusive MC events. The main 
backgrounds come from the decays of "0(3686) — > p^K^K* + p°K* ± K* T — > TT ± ir T ir°KgK ± ir T , which don't form a 
peak in the M K o sK +^- spectrum. The background from the e + er — > qq continuum process is studied by using data 
collected at the center-of-mass energy of 3.65 GeV. Continuum backgrounds are found to be small and uniformly 
distributed in the mass spectrum. 

An unbinned maximum-likelihood fit to M K o K + 7r ~ is applied to determine the number of wA(1440) and a; /i (1285) 
events. The fit includes three components: A(1440), Ji(1285), and background. Both A(1440), and /i(1285) are 
represented by Brcit-Wigner (BW) functions convoluted with a Novosibirsk mass-resolution function and multiplied 
by an efficiency curve. Both the mass resolution and the efficiency curve arc determined from MC simulations. The 
background shape is described by a 2nd-order Chebychev polynomial function. The mass and width of /i(1285) are 
fixed at the known values [14, 27], while those for the A(1440) are allowed to float in the fit. The solid curve in 
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M„.„„ (GeV/c 2 ) M K „ KV (GeV/c 2 ) 

FIG. 1: (a) Distribution of M^+ n - versus M^+^-^o, where the boxes represent the K s and w signal region and sideband 
regions, (b) The K%K + ir~ invariant -mass distribution for ^>(3686) — > LjKgK + iv candidate events. The shaded histogram is 
the background contribution estimated from the Kg and ui sidebands minus that of the corner ranges. 

Fig. 1(b) shows the fit results. The goodness-of-fit is % 2 /ndf = 81.0/75 = 1.08, which indicates a reasonable fit. The 
fit yields 109.2 ± 18.0 events for the X (1440) with a statistical significance of 9.5a, and 21.6 ± 7.0 events for the 
/i(1285) with the statistical significance of 4.5er. The significance is determined by the change of the log-likelihood 
value and the degrees of freedom in the fit with and without assuming the presence of a signal. The mass and width 
of the X(1440) are determined to be M = 1452.2 ± 5.2 MeV/c 2 and V = 51.6 ± 12.1 MeV/c 2 . The contribution of 
background estimated from the sideband is 13.4 ± 6.1 for the X(1AAQ) and 0.1 ± 2.4 for the /i(1285). In the analysis, 
these backgrounds arc ignored. 

With the MC-determined efficiencies (where the final-state particles distribute uniformly over the phase space) of 
(10.41 ± 0.14)% for ^(3686) wX(1440) and (10.88 ± 0.15)% for ^(3686) -> w/i(1285), the product of branching 
fractions are calculated to be 

23(^(3686) -> wX(1440)) ■ B(X(1440) -» K° s K + t:-) = (1.60 ± 0.27 (stat.)) x 10~ 5 , 

S(^(3686) —¥ w/i(1285)) • B(/ x (1285) -> A'g/v+O = (3.02 ± 0.98 (stat.)) x 10~ 6 . 
Systematic uncertainties are discussed in Section VI. Since the /i(1285) peak is not significant, an upper limit at a 
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90% confidence level (C.L.) on N sig is determined using a Bayesian method [14] by finding the value N^g such that 

mUP 

f n 319 LdN siq 

^— S J± = 0.90, 

Jo LdN sig 

where N sig is the number of signal events, and L is the value of the likelihood as a function of N sig . An upper limit 
at the 90% C.L of 31 /i(1285) events is obtained. 

B. Branching fractions of ^(3686) -> ujK* + K~ + c.c. and ^(3686) -s- ujK*°K° 



Figure 2 shows the M K o n ± and M K ^ 7T ± distributions, where the A*(892) and a peak at 1.43 GeV/c 2 are evident 
in both distributions. The solid lines in Fig. 2 show the results of an unbinned maximum likelihood fit with three 
components: A*(892), A2(1430), and background. Both the A*(892) and ^(1430) are described by acceptance- 
corrected BW functions. The BW function [28-30] is 

MT(s) 

*BW(S)- (s 2_ M 2)2 +M 2 r(s) 2' W 

where T(s) = r(^-) 2 (^) 2i+1 , M and T are the K* mass and width, q is the K° s (A=F) momentum in the K* rest- 
frame, go is the q value at s — M, and L is the relative orbital angular momentum of Kgir^ (K^tt^ 1 ). Here the 
A*(892) and A"2(1430) peaks are described with P-wavc (L = 1) and D-w&ve (L = 2) BW functions, respectively. 
The mass and width of the A* (892) are floating, and those of A2(1430) are fixed to the world average values [14]. 
The backgrounds are described by the function [31, 32] 

Fbg(s) = (s - m t ) c e- ds - es2 , (2) 

where nit is the threshold mass for K^tt^ (K^tt^) and c, d and e are free parameters. 

The fit to the M K o^± distribution yields 502.0 ± 56.4 A* ± (892) events, and 128.5 ± 30.0 A'* ± (1430) events with a 
statistical significance of 4.4<r. The mass and width of the A* ± (892) determined in the fit are 888.0 ±2.5 MeV/c 2 and 
48.0 ± 6.5 MeV/c 2 , respectively. The fit to M KTv ± yields 446.2 ± 47.4 A*°(892) events, and 164.2 ± 34.2 A*°(1430) 
events with a statistical significance of 4.6<7. The mass and width of the A*°(892) are determined to be M = 893.9±2.2 
MeV/c 2 and T = 45.0 ± 5.7 MeV/c 2 . When the peak at 1430 MeV/c 2 is fitted to an S-wave A*(1430) or a P-wave 
A* (1410), the fit qualities degrade for both fits. The A 2 *(1430) can be distinguished from A o *(1430) and A'* (1410) 
with log-likelihood values worse by 173.9 and 173.7, respectively, in the A^7r ± decay, while the log-likelihoods are 
worse by 170.6 and 169.7, respectively, in the K^n^ decay. 
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The peaking backgrounds for ^sT* (892) and ^(1430) are studied using the Kg and u sidebands. In K^tt^, the 
contributions from background events are estimated to be 105.6 ± 21.6 for if* ± (892), which are subtracted from 
the total signal yield, while no evident peaking backgrounds for K^il^SO) are seen. Combining the numbers of 
signal events with the detection efficiency of (9.58 ± 0.08)% and (9.18 ± 0.07)% determined from MC simulation of 
-0(3686) — > ojK* + (892)K~ and 0(3686) — > ujK^ (1430) JsT - , respectively, their corresponding branching fractions are 
determined to be 

i?(0>(3686) loK*+(892)K-) = (1.89 ±0.29 (stat.)) x 10~ 4 , 
S(V>(3686) ->■ u)Kl + {lAZQi)K-) = (6.39 ±1.50 (stat.)) x 10~ 5 . 




M k°, ( QeV/c2 ) M K3I (GeV/c 2 ) 

FIG. 2: (a) The M K o 7I ± distribution, (b) The M XTlr ± distribution. Solid curves show the fit results; dashed lines are for 
signals; and dotted lines are for backgrounds. 

In the K^ir* mo dc, after subtracting 90.2 ± 18.3 background events for K*°(892) and 48.5 ± 23.2 for A' 2 *°(1430), 
and using the MC-determined detection efficiencies of (9.66 ± 0.08)% and (9.08 ± 0.07)%, respectively, the branching 
fractions are 

23(0(3686) -> ujK*°(892)K°) = (1.68 ±0.25 (stat.)) x 10~ 4 , 



6(0(3686) -> ojKZ°(U30)K ) = (5.82 ± 2.08 (stat.)) x 10" 5 . 
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V. 0(3686) —¥ ujK + K~ tv° 

Candidate events must have four charged tracks with zero net charge, and four or more photons. At least two 
charged tracks should be identified as kaons. When more than two charged tracks are identified as kaons the two 
oppositely charged tracks with the largest Prob(-ftT) are chosen. To reject the two backgrounds with 37 or 57 in the 
final state, the xlc °f a four-constraint kinematic fit in the hypothesis of ^(3686) — > K + K^tt + tt^A^ is required 
to be less than those of the K + K~ir + ir~3-f, and K + K~Tr + ir~5-f hypotheses. For all possible combinations, a six- 
constraint kinematic (6C) fit (besides the initial 0(3686) four-momentum, two ir° masses are also used as constraints) 
is performed, and the one with the least \\c ^ s chosen, and x\c < 100 i s required. The combination of 7r + 7r~7r 
with the invariant mass closest to the nominal w mass is selected as an lu candidate; the invariant mass M^+^-^o 
must be in the region 0.743 GeV/c 2 < M w +„-„o < 0.823 GcV/c 2 . Background events from 0(3686) — > tt + it~ J/ip and 
0(3686) 7r 7r° J/0 decays are removed by requiring \M™ C ° U - m J/v ,| > 0.007 GeV/c 2 and |Af;^° li - m J/4 ,\ > 0.06 
GeV/c 2 . Background events from 0(3686) — > 77J/0 are rejected by requiring \M^ cml — mj/^\ > 0.05 GeV/c 2 . 

A. Branching fractions for 0(3686) -)■ uX(lUO) -)■ uiK + K~n° 

After the above event selection, the distribution of Mk+k-^o is shown in Fig. 3(a), where the A(1440) is clearly seen 
and the /i(1285) is evident. Non-w background, estimated from the to sideband (0.06 GeV/c 2 < |A/„.+„— ^0 — m w < 
0.10 GeV/c 2 ), is shown as the shaded histogram and doesn't form a peak in the M^+k-tt spectrum. 

The results of an unbinned maximum-likelihood fit to M K + K -^o, which is similar to the M K o K7I fit described in 
Section IVA, is shown as a solid line in Fig. 3. The goodness-of-fit is x 2 /ndf = 63.6/67 = 0.95, which indicates 
a reasonable fit. The fit yields 81.8 ± 14.7 A(1440) events with a statistical significance of 9.3a, and 9.5 ± 5.3 
/i(1285) events with a statistical significance of 3.2cr. The mass and width of the X(1440) determined from the fit 
are M = 1452.7 ± 3.8 MeV/c 2 and V = 36.8 ± 10.5 MeV/c 2 , respectively. The non-w contributions estimated from 
the Lo sideband are 2.8 ± 3.9 events for the A(1440), and 0.1 ± 1.5 events for /i(1285), which are neglected in the 
branching fraction measurements. 

With the detection efficiencies determined from a phase space distributed MC simulation of ■0(3686) — > cjA(1440) 
and -0(3686) cj/i(1285), which are (7.92 ± 0.13)% and (8.02 ± 0.13)% respectively, the products of branching 
fractions are determined to be 
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60(3686) -)• wX(1440)) ■ S(X(1440) -> K+R-n ) = (1.09 ± 0.20 (stat.)) x 10~ 5 , 

#0(3686) -4 w/i(1285)) • B(/ x (1285) K+R-n ) = (1.25 ± 0.70 (stat.)) x 10~ 6 . 

Using the Bayesian approach, the upper limit on the branching fraction for 0(3686) — > cj.fi (1285) is determined. 
The upper limit on the number of /i(1285) events is 15 at the 90% C.L.. 

B. Branching fractions for i/>(3686) — s- u)K* + K~ + c.c. 

The M K ±^o distribution (both M K + n o and M K - n o are included) is shown in Fig. 3(b), where the if* (892) and 
A",*(1430) are clear. Using the same functions as described in Section IVB, an unbinned maximum likelihood fit 
to M K±7r o is performed. The fit yields 678.8 ± 65.3 ijf* ± (892) and 142.8 ± 39.0 K^(U30) events with a statistical 
significance of 4.5cr. The mass and width of A'* ± (892) are determined to be M = 889.6±2.1 MeV/c 2 and V = 49.2±5.5 
McV/c 2 . 

The background contributions from non-cj processes estimated by fitting M K ±„o for events in the u sideband 
are (144.2 ± 24.8) events for A'* ± (892), which is subtracted from the total K* ± (892) yield, while no evident peak 
contributions for if ^ (1430) are present. 

The efficiencies determined from MC simulation are (7.48±0.07)% and (7.70±0.07)% for 0(3686) ->• loK* + {892)K~ 
and '0(3686) — > wif2 + (1430)if~, respectively. The branching fractions are determined to be 

80(3686) ->• luK* + (892)K ~) = (2.26 ± 0.30 (stat.)) x 10" 4 , 
£(0(3686) -»• ujK* + (U30)K-) = (5.86 ± 1.61 (stat.)) x 10" 5 . 

VI. SYSTEMATIC UNCERTAINTIES 

Systematic errors in the branching fraction measurements mainly come from the number of ■0(3686) events, tracking, 
particle identification, photon reconstruction, K s reconstruction, kinematic fit, background estimation, signal shape 
and detection efficiency. 
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FIG. 3: (a) The M K + 

K- n ° invariant-mass distribution for ^(3686) — > u>K + K tt° candidate events. The shaded histogram 
is for non-o; background events estimated from the u sideband, (b) The combined mass spectra of K + -k° and K~tt° for 
0(3686) — > ljK* + K~ candidate events. Solid curves are the fitting results; dashed lines are for the signal; and dotted lines are 
for the background. 

The uncertainty on the number of 0(3686) is 0.81% as reported in Rcf. [35]. The tracking efficiencies have been 
checked with clean channels such as J/ip — > pn, J/ip — > Tr + 7r~pp, 0(3686) — > 7r + 7r~ J/0 (J/ip —> l + l~) and J/ip — > 
A*°(892)Ag (K*° -> A^tt^, K s -> tt+tt"). It is found that the MC simulation agrees with data within 1% for 
each charged track [33]. A 4% systematic error due to the MDC tracking efficiency is assigned for both -0(3686) — > 
ujKgK + ir~ and 0(3686) — > ujK + K~ir° . The uncertainty associated with the kaon identification has been studied 
with a clean kaon sample selected from J/ip —> K s K + ir~ [33]. The data-MC efficiency difference, which is 2% per 
track, is assigned as the systematic error for the kaon identification. 

The uncertainty in the K s reconstruction efficiency [34] , which includes the geometric acceptance, a pair of pion 
tracking cfficcncies and the K s secondary vertex fit, is estimated with the decay of J/0 — > K*°K°, and the MC-data 
difference of 3.5% is taken as the systematic error for the K s reconstruction in the decay 0(3686) — > ujK s K + it~ . The 
photon detection efficiency is studied with J/ip —> p°7r° , and the difference between data and MC simulation is about 
1% [33] per photon. A 2% (4%) systematic error for photon efficiency is used for the decay -0(3686) — > ujKgK + ir~ 
(0(3686) -> ojK+R-t: ). 

The uncertainties associated with the kinematic fit arc estimated by using 0(3686) — > tt + tt~J/tJj with J/ip — > 
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KgK + ir~Tr° or J/0 —> K + K~tt tt° events. The efficiencies are obtained by comparing the number of signal events 
with and without the kinematic fit performed for data and MC simulation separately. A data-MC difference of 5.4% 
is found in J/ip KgK + ir~ 7r° and 3.2% in J/0 — > K + K~ir Q n° . The differences are taken as the systematic errors. 

The uncertainty in the 6(0(3686) — > a;A(1440)) measurement due to the background shape is estimated by varying 
the background function from a 2nd-order polynomial to a 3rd-order polynomial in the fit to M K j^^. The changes 
in X(1440) and /i(1285) yields, are 5.1% (3.2%) and 2.5% (2.2%), respectively, in the decay 0(3686) ujK° s K+tt- 
(-0(3686) — > luK + K~tt°). The uncertainty due to the fit range is estimated by repeating the fits in the range [1.2,2.4] 
GcV/c 2 , and the differences of 6.1% (3.3%) and 3.3% (6.0%) in the branching fractions for X(1440) and /i(1285) in 
the decay -0(3686) — > u>KgK + Tr~ (-0(3686) — > uiK + K~tt°) are assigned as systematic errors. 

The uncertainty in the 23(0(3686) — > o;X(1440)) measurement due to the signal shape is considered to come from 
the mass resolution, mass shift and phase space factor. By varying the mass resolution by ±0.5 MeV/c 2 from the MC 
expectation, the differences in A(1440) yield, 1.0% for KgK + ir~ and 1.1% for K + K~ir°, are assigned as systematic 
errors respectively; the difference in the /i(1285) yield can be ignored. By floating the mass of /i(1285), the changes, 
1.0% for KgK + ir~ and 2.1% for K + K~ir°, are taken as the systematic errors. After taking the phase space factor into 
account, which depends on the momentum of the X(1440) in the 0(3686) rest frame and the relative orbital angular 
momentum between the w and A(1440), the differences in the fitting results of A(1440) and /i(1285) are found to 
be 2.7% and 1.0% in the decay 0(3686) ->■ ujK%K + -k- , and 0.5% and 2.1% in the decay 0(3686) ->■ uK+R-ir , 
respectively, which are taken as the uncertainties. 

The selection efficiencies are determined from phase space distributed MC simulations of the 0(3686) decay. The 
uncertainties in the selection efficiencies are estimated by using efficiencies obtained from MC samples that include 
intermediate states, or the angular distribution associated with the A(1440). There is a small difference (2.5% for 
KgK + ir~ and 3.2% for A" + A"~7r°) from the A*A intermediate process. The detection efficiencies are also checked 
by generating 0(3686) — > wA(1440) MC events assuming X(1440) is a pseudoscalar meson, and the differences for 
A(1440) -> A'gA+TT" and A(1440) -> A+ A^tt are 10.0% and 11.3%, respectively. To be conservative, the differences 
of 10.0% and 11.3% are taken as the systematic errors for the A(1440) — » K%K + %- and A(1440) A+A^tt , 
respectively. 

The uncertainties in the £>(0(3686) — > uiK*K) measurement due to the fit range are estimated to be 7.6% (8.9%), 
3.6% (9.1%) and 3.6% (11.6%) for the decay A* ± (892)(A* ± (1430)) —¥ K^, A"*°(892)(A 2 *°(1430)) K^n ± and 
AT* ± (892) (K^ (1430)) -> A^tt , respectively. The differences by changing the sideband range (0.014 GcV/c 2 < 
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TABLE I: The systematic errors (%) of 23(^(3686) ->■ loX -)■ w^Atf). 



Sources 


LjKgK TV 


UjK K TV 




A (1440) /i(1285) 


A (1440) /i (1285) 


The number -0(3686) events 


0.8 


0.8 


MDC tracking 


4 


4 


Particle identification 


2 


4 


Kg reconstruction 


3.5 




Photon efficiency 


2 


4 


Intermediate decays 


0.8 


0.8 


Kinematic fit 


5.4 


3.2 


Background uncertainty 


6.1 3.3 


3.3 6.0 


Signal shape 


2.9 1.4 


1.2 3.0 


MC Eff. Uncertainty 


10.0 


11.3 


Total 


14.6 8.9 


14.2 10.3 



|M w+ff - - m K o \ < 0.022 GcV/c 2 or 0.08 GcV/c 2 < \M„+ n -^o - m u \ < 0.12 GcV/c 2 ) are estimated for the above 
decay modes. For if* (892), the differences in the fitting results with or without the presence of if|(1430) or by 
replacing i^2(1430) with ifg(1430) (or i4T*(1410)) are also estimated. The above largest differences are taken as the 
systematic errors of the A"* (892). 

All the contributions are summarized in Table 1 and Table II. The total systematic uncertainty is given by the 
quadratic sum of the individual errors, assuming all sources to be independent. 

VII. DISCUSSION 

As the X(1440) and /i(1285) are observed in both K^K^ir^ and K + K~tt° final states, a simultaneous maximum- 
likelihood fit is performed to the mass spectra to extract a more precise determination of the resonant parameters 
and branching fractions. The fit includes three components, the X(1440), /i(1285), and background, as used in the 
fit to each individual mode in Section IV A and Section V A. The fit, shown in Fig. 4, has a x 2 /ndf = 72.9/70 = 1.04, 
and the statistical significances of the AT(1440) and /i(1285) are 13.3tr and 5.4a, respectively. The mass and width of 
the X(1440) from the fit are M = 1452.7 ± 3.3 McV/c 2 and T = 45.9 ± 8.2 McV/c 2 . The yields of the AT(1440) events 
arc 111.4 ± 17.2 in the K^^tt^ mode and 82.4 ± 13.5 in the K+R-tt mode, while those of the A (1285) events arc 
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TABLE II: The systematic errors (%) of B(i/>(3686) -> luK*K). 



Sources 




LjKgK TV 




ujK + K 




TV 




K (892) K 2 


(1430) K (892) K 2 


(1430) 


A (892) K 2 


±(1430) 


The number of t/>(3686) events 




0.8 




0.8 




MDC tracking 




4 




4 




Particle identification 




2 




4 




Kg reconstruction 




.3.5 








Photon efficiency 




2 




4 




Intermediate decays 




0.8 




0.8 




Kinematic fit 




5.4 




3.2 




Background uncertainty 


7.6 


8.9 3.6 


9.1 


7.2 


11.6 


Total 


11.2 


12.1 9.0 


12.3 


10.6 


14.0 



23.1 ± 7.1 in the K^K^-k^ mode and 8.7 ± 4.6 in the K + K ~n° mode, in good agreement with the separate fits to 
the two modes, as shown in Table III. Combining the observed numbers of signal events with efficiencies and taking 
properly into account the Clebsch-Gordan coefficients, the branching fractions are 

£(V>(3686) -> ujX(IUO)) ■ B(X(UA0) ->■ KKn) = (5.48 ± 0.61 (stat.) ± 0.86 (sys.)) x 10" 5 , 

and 

B(^(3686) —¥ w/i(1285)) • S(/i(1285) -> KK-k) = (8.78 ± 2.33 (stat.) ± 0.96 (sys.)) x 10" 6 . 

The first errors are statistical and the second ones systematic. In calculating the systematic errors, the correlations 
between the errors in the two modes shown in Table I are properly taken into account. 

To examine the spin-parity of the events at around 1.4 GeV/c 2 , we try to measure the polar angle distribution 
(denoted as cos^x) of the X(1440) in the '0(3686) rest frame. The | cos9x \ distribution is divided into five bins in the 
region of [0, 1]. In each bin, the combined KK-k mass spectrum of the two decay modes is fitted. For X(1440), the 
mass and width are fixed to those of the combined mass spectrum fitting for the whole angular range; for /i(1285), 
the mass and width are fixed at the known values. By repetition of the mass fit in five bins of | cos#x|, the number 
of X(1440) events can be obtained. Figure 5 shows the polar angular distribution for signal yields, where the errors 
are statistical only. 

The angular distribution is fitted to 1 + acos 2 #x, as shown in the solid line in Fig. 5, and a = 0.58 ± 0.64 is 
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FIG. 4: The simultaneous maximum-likelihood fit to the mass spectra for KgK tc^ (a) and K + K ir° (b). Points with error 
bars are data, and the curves show the best fits. 

obtained with a probability of 29%. Around 1.44 GeV/c 2 , there are two known resonances, namely the 77(1440) with 
J PC = CT+ and the / x (1420) with J PC = 1++ [14]. The present statistics arc not sufficient to establish a 1 + cos 2 Ox 
behavior as a pseudoscalar meson candidate of ?y(1440), but the data clearly favor a = 1 or a = over a = —1, as 
can be seen from Fig. 5. 



VIII. SUMMARY 



With a sample of 106 million -0(3686) events, the decay of -0(3686) — > ujKKit is studied for the first time. In 
addition to the mass enhancement (X(1440)) around 1.44 GeV/c 2 , the /i(1285) is observed in the mass spectrum of 
KKtt. From investigating the X(1440) polar angle distribution, the present statistics are not sufficient to establish 
a 1 + cos 2 Ox behavior as a pseudoscalar meson candidate of r/(1440), but the data favor a = 1 or a = over 
a = — 1. The product branching fraction upper limit of /i(1285) in each individual mode is also presented after 
taking into account the systematic error by dividing by a factor (1 — 5 sys ), where the S sys is the systematic error for 
the corresponding decay. Also the branching fractions of 0(3686) — > u>K*K for the charged and neutral mode are 
measured for the first time. The observed if*(1430) favors ^(1430) over ^(1430) and If* (1410). The numbers of 
observed events, detection efficiencies and branching fractions (or upper limits) are summarized in Tables III and IV. 

To compare with the 12% rule, Tables III and IV also include the corresponding J/0 branching fractions [36], as 
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X 

<x> 




|cosO x | 

FIG. 5: The | cos Qx | distribution with the fit results; solid line is for 1 + a cos 2 Qx with a — 0.58 ± 0.64; dash-dotted line is for 
1 + cos 2 Ox; dotted line is for 1 + • cos 2 8x; and dashed line is for 1 — cos 2 Qx- The expectation distributions are multiplied 
by the efficiency. 

TABLE III: The branching fractions and upper limits (90% C.L.) for t/j(3686) — > uX — > uuKKtt decays. Results for correspond- 
ing J/ip decays [36] and the ratio Qh = jy^ 6 ^)^ are a ^ so gi ven - 



decay mode 


N St g 


e(%) 


B(V>(3686))(xlO- 5 ) 


B(J/V)(xl0- 5 ) 


Q h (%) 


u;X(1440) -> 




109 ± 18 


10.41 ±0.14 


1.60 ±0.27 ±0.24 


48.6 ±6.9 ±8.1 


3.3 ± 1.1 




LOK+K-TT 


82 ± 15 


7.92 ±0.13 


1.09 ±0.20 ±0.16 


19.2 ±5.7 ±3.8 


5.7 ±2.5 


— » 


ljKKiy 






5.48 ±0.61 ±0.86 






u) fi (1285) - 


ujK%K + ti- 


21.6 ±7.0 


10.88 ±0.15 


0.302 ± 0.098 ±0.027 










< 31 


10.88 ±0.15 


< 0.478 








LOK+R-TT 


9.5 ±5.3 


8.02 ±0.13 


0.125 ±0.070 ±0.013 










< 15 


8.02 ±0.13 


< 0.221 






— s 


ljKKty 






0.878 ± 0.233 ± 0.096 







well as the ratio Qh- The data show that -0(3686) decaying into wX(1440) and luK*(892)K are suppressed by a factor 
of 2 - 4. 
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TABLE IV: The branching fractions of ^(3686) -s- uK*K. 



decay mode 


Y 3i9 


e(%) 


B(^(3686))(xl0 -5 ) 


S(J/V0(xlO 5 ) 


Qh{%) 


ujK*+(892)K- - 


uK%K + ir~ 


396.4 ± 60.4 9.58 ±0.08 


18.9 ± 2.9 ± 2.2 


310 ± 34 ± 56 


6.1 ± 1.8 




■> cjK+K~ty° 


534.6 ± 69.9 


7.48 ± 0.07 


22.6 ± 3.0 ± 2.4 


327 ± 51 ± 68 


7.0 ±2.2 


uK^ + (1430)K~ 


-> uiKgK + n~ 


128.5 ±30.0 


9.18 ±0.07 


6.39 ±1.50 ±0.78 








-> ljK+K+tt 


142.8 ± 39.0 


7.70 ± 0.07 


5.86 ±1.61 ±0.83 






luK*°(892)K° - 


► ujK%K + -k~ 


356.0 ± 50.8 


9.66 ± 0.08 


16.8 ± 2.5 ± 1.6 


310 ± 34 ± 56 


5.4 ± 1.5 


uKf(U30)K° - 




115.7 ±41.3 


9.08 ±0.07 


5.82 ±2.08 ±0.72 
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